Технические науки 


УДК 625.032.3:681.3.06+06 


ДИНАМИЧЕСКИЙ МОНИТОРИНГ ТРИБОТЕРМОДИНАМИКИ 
ФРИКЦИОННЫХ МОБИЛЬНЫХ СИСТЕМ 


А.Л. ОЗЯБКИН 
(Ростовский государственный университет путей сообщения) 


Выдвинута и доказана на практическом примере гипотеза экспериментальной триботермодинамики, позво- 
ляющая устанавливать значение температуры вспышки на фактической площади касания. 

Ключевые слова: Ффизико-математическое моделирование, трибоспектральная оценка, температура 
вспышки, схватывание, поверхностная и объемная температуры. 


Введение. Наиболее напряженные условия взаимодействия подсистемы «колесо — рельс» скла- 
дываются на сортировочных горках и при вписывании товарного подвижного состава в криволи- 
нейные участки пути малого радиуса с отличным от нуля непогашенным ускорением, обусловлен- 
ным пониженной его скоростью движения по сравнению с нормативной. В результате такого дви- 
жения подвижного состава под действием потенциальных сил гребни колесных пар долговремен- 
но контактируют с боковой гранью головки внутреннего рельса. Значительные нагрузки в сочета- 
нии со значительными скоростями скольжения приводят к повышению температуры в контактной 
области колеса с рельсом (ее значение может быть сопоставимо со значением температуры плав- 
ления материала колеса), что способствует изменению физико-механических свойств фрикцион- 
ных пар и неблагоприятно сказывается на стабильности динамических характеристик. Возникаю- 
щие при движении железнодорожного вагона динамические силы являются следствием колеба- 
тельных и других видов неравномерного движения инерционных масс, составляющих рассматри- 
ваемую механическую систему [1] и соединенных между собой упругодиссипативными связями. 
Математическое моделирование такой системы [2] всегда осуществляется с учетом ряда допуще- 
ний, что может приводить к неадекватным результатам моделирования, особенно, если это каса- 
ется триботермодинамики. 
Моделирование триботермодинамики с учетом допущений. В качестве примера реализа- 
ции динамического мониторинга рассмотрим грузовой вагон и его подсистему «колесо — рельс» 
для идентификации его триботермодинамики. Выполним расчет максимальной температуры в з0- 
не трения колеса грузового вагона с рельсом при роспуске вагонов с сортировочной горки с баш- 
мачным торможением по известной методике [3]: 

Ох = Оу + Ош» + Ове, (1) 
где 9„.х — суммарная температура в контакте, °С; 9’ — объемная температура, °С; © - поверх- 

ностная температура, °С; ®9„„ - температура вспышки, °С. 
Объемная температура определяется по формуле: 


_ @-о.) (2) 
5’усЬЕ. 

где о. — коэффициент распределения тепловых потоков между элементами пары трения; 
Й’,- полная работа трения за время торможения, Дж; 5 — номинальная площадь контак- 
та, м”; у — плотность материала (7800 кг/м“); с — удельная теплоемкость (503 Дж/(кг.град)); 

р — толщина рассматриваемого элемента пары трения, м; #.— время торможения, с. 
Для расчета номинальной площади контакта воспользуемся формулой для случая контак- 
та двух цилиндров с параллельными осями. Допустим, что площадь контакта имеет форму прямо- 
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угольника со сторонами /хБ. Здесь [ — ширина горизонтального участка рельса (49 мм), Б — длина, 


определяемая как 
У и 
ЕЮ В 


где Р — статическая нагрузка вагона на рельс (133,5 кН); Е — модуль Юнга (2,1.10*' Па); А, - ра- 
диус колеса вагона (0,475 м); К› — радиус рельса (‹о, м). 

Тогда Р=3,78 мм, а номинальная площадь контакта 5=185,14 мм”. Определим полную ра- 
боту сил трения за время торможения 
И’ = Е-А[ = 24965.12 = 299580 Дж, 

где А/ — расстояние между противоползунными конструкциями ППС-12Д, предназначенными для 
предотвращения образования ползунов колесных пар; Р-— равнодействующая сила, дейст- 
вующая на колесо вагона при движении с сортировочной горки: 

Е = те зш о — пе с0$ 4. 

Здесь т — распределенная масса вагона, приходящаяся на одно колесо (13,35 т); о — угол накло- 
на горки по отношению к горизонту; Г-— коэффициент трения (0,24); = - ускорение свобод- 
ного падения. 

Определим коэффициент распределения тепловых потоков: 


1 1 
а 0 
[4-22 2 1 
с, | 0,18 
где ат, а> —- коэффициенты температуропроводности колеса и рельса (так как физические свойст- 
ва колеса и рельса приблизительно равны, то а! = а›=17,65:10® град-м”/с); Ь, — размер 
пары трения в направлении теплового потока (высота рельса 180 мм); Б› — диаметр колеса 
(950 мм); сти с> — удельная теплоемкость колеса и рельса. 
После нахождения всех параметров, объемная температура [см. формулу (2)] составит: 
9 — 0,16.299580 
’ 185,14.10°.7800-503.0,95.4,3 
Определим значение средней поверхностной температуры: 


Е } 


"м + Рот 








= 16, [5 °С. 


ПОВ у [. 5 


где ^, — коэффициент теплопроводности (57,7 Вт/(м.-град)); 5-— номинальная площадь поверхности 
трения, м*; ху и ти — безразмерные параметры работы и мощности трения, изменяющиеся 


(3) 


по закону торможения вагона т, = (0,5%). [3 В (0,5%), т — безразмерный коэффициент 


времени торможения, т = {/ д; Е — время трения; Ео — число Фурье, значение которого за- 
висит от эффективной и действительной толщины элементов пары трения: если Б <р,, то 


[ТА 
Ро те, а если Б>Ь. , то Ео = 1/3; Б., —- эффективная толщина материала, участвующего в 


теплопоглощении, которую можно вычислить согласно выражению: 
Ь. =1,73./аё, =1,73\/17,65.10°.4,3 = 0,015 м. 


После нахождения всех параметров, выражение (3) для поверхностной температуры коле- 
са при трении примет вид: 
0,16.299580. 0,015 


м 57,7-4.3.185.14- 10-3 (ити (ити). 
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Результаты расчетов сведены в таблицу. 
Расчетные величины, определяющие поверхностную температуру 


Л СИИ 
р ро в 
р р ры 
ПС ЗВ МЮ И С М И И 
ПВ ИИ И ИИ М. С ИЯ 
ОС ПИ 

о 


р р рр 


Получено, что поверхностная температура уже тогда достигает высоких значений, т.е. 
температуры плавления материала, когда большая часть активных микрообъемов на фактической 
площади касания подвержена не упругим (по закону Гука), а пластическим деформациям. В дан- 
ном случае расчет температуры вспышки не приводится, так как он даст еще большие значения. 
По нашему мнению, в выражении (1) должны указываться не абсолютные значения температур, а 
их приращения как функции градиента температуры. Кроме того, указанный способ мониторинга 
триботермодинамики применять невозможно, так как он не учитывает взаимосвязь и взаимовлия- 
ние динамической системы, подсистемы фрикционного контакта и условий окружающей среды. 
Гипотеза экспериментальной триботермодинамики. Нами предложен способ уточнения 
выражения (1), позволяющий проверить экспериментальную модель триботермодинамики фрик- 
ционной мобильной системы, в частности железнодорожного вагона, с применением метода фи- 
зико-математического моделирования [4]. Процессы в объекте исследования и модели при дина- 
мическом подобии должны быть тождественны и описываться одинаковыми дифференциальными 
уравнениями [4, 5]. В качестве условия динамического подобия механической подсистемы вагона 
принимаем равенство частот колебаний модели и объекта исследования. 

Основным критерием достоверности модельного эксперимента следует считать реализа- 
цию в условиях модели одинакового с натурным образцом вида изнашивания поверхностных 
слоев и интенсивности износа контактирующих тел. Для реализации условий динамического по- 
добия процессов на поверхности фрикционного контакта при трении необходимо выполнить сле- 
дующее [4]: 





в ме ое! 


Сс ес. =С=С С, Сс С, 
где масштабы подобия: С», С, — собственных и затухающих частот колебаний соответственно; 
Ст- периода колебаний; С, — высоты микрообъема; С. - радиуса закругления микрообъема; 
С. - амплитуды деформации микрообъема; Су„ — скорости относительного скольжения; 
С, - контактного давления; Су — фактической площади касания; С+ — силы трения; 
Суи- нормальной нагрузки; С» — массы микрообъема; Сс — жесткости микрообъема; Св - эк- 


вивалентного коэффициента демпфирования; С) — геометрических размеров; С’, — момента 
инерции микрообъема. 
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Фрикционный узел трения «колесо — рельс» является составной частью фрикционной мо- 
бильной системы «подвижной состав -— путь». Масса колеса т, и приведенная масса рельса то со- 
вершают сложное взаимное перемещение и являются составными частями механической системы 
(рис.1). 

Теоретический расчет или экспериментальное исследование прямыми методами теплового 
состояния А-Б и Б-В области на сегодняшний день являются проблемными, так как определяются 
текущими значениями шероховатостей контактирующих поверхностей при их взаимном смеще- 
нии, значением фактической площади касания, которые меняются на порядок и более при отно- 
сительном скольжении трущихся тел [6] и определяют триботехнические характеристики фрик- 
ционных систем. Значения объемных температур в В-Г области, т.е. ее градиентное изменение, 
можно определить на базе существующих математических моделей. 
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Рис.1. Модельное представление фрикционной системы «колесо -— рельс»: Р’у- нагрузка; [’, — сила 
трения; / — коэффициент трения; 7711 — масса колеса; 71> — приведенная масса рельса; © - темпера- 
тура; 9»ст — максимальная температура, реализуемая на наибольшей высоте А,; ту — объемная 
масса; 77, — поверхностная масса; тт„ — масса контактирующих микронеровностей поверхности; 
Ю. — высота неровности профиля по десяти точкам; А, — наибольшая высота профиля; Эном — но- 
минальная площадь поверхности; фк — фактическая площадь касания поверхностей трения; 


'/ - плотность материала; Р - усилие; 2 — сближение поверхностей трения; .4., 4. — перемещение 
колеса относительно рельса в горизонтальном и вертикальном направлениях соответственно 
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Выдвигается гипотеза, что максимальная температура в контакте ©... [см. формулу (1)] не 
складывается из трех отдельных температур (объемной, поверхностной и температуры вспышки), 
как указывается в [3], а является одной температурой, имеющей три зоны состояния [7]: 

— область А-Б асимптотического изменения объемной температуры от максимальных зна- 
чений (температура вспышки) до максимальной поверхностной, значения данной температуры 
могут достигать и превышать температуру плавления контактирующих материалов; 

— область Б-В перехода объемной температуры от максимальной поверхностной к макси- 
мальной объемной в сплошном теле; 

— область В-Г асимптотического перехода от максимальной температуры в сплошном теле 
до ее низших значений. 

В качестве реализации гипотезы триботермодинамики принимается постоянная величина 
произведений масс активных микрообъемов на показатель асимптоты градиента объемной темпе- 
ратуры, т.е. 

т, =т “> = Тб. (4) 

Для подтверждения выдвинутой гипотезы, увеличения достоверности и адекватности экс- 
периментального определения значений физических параметров, установления законов статисти- 
ческой повторяемости результатов, имеющих вероятностную природу, проведения математиче- 
ского планирования эксперимента, установления рациональных значений параметров и переноса 
результатов испытаний на реальный объект исследования методом анализа размерностей выво- 
дятся константы подобия физических величин, существенно влияющих на триботермодинамику 
трибосистемы «колесо — рельс». Функциональная зависимость поверхностной температуры ©; от 
физических параметров имеет вид критериального уравнения подобия [5]: 


М а ЛОГ Гб Е т’ ВИ РГР И 

5 Ум’ е’ РМ мм М °М’, 
5-7 ГГ 52 ТР ог 20 у ОР АЙ 
МУ 125,’ МУ’ МУ °М’М’ М Г 
где базисные факторы: ® — максимальная температура на фактической площади касания, °С; 
№ — нормальная нагрузка колеса железнодорожного вагона на рельс, Н; Г -— скорость каче- 
ния, м/с; /[— геометрический масштаб, м; граничные условия; К — частота собственных коле- 
баний, с '; д — контактное давление, Па; Л® - градиент температуры, К/м; прочие факторы: 
1 — время трения, с; Г, -путь трения, м; 5 — фактическая площадь касания, м”; Е’ - сила тре- 
ния, Н; т — масса колеса, кг; В — коэффициент демпфирования, Н.с/м; р — удельный вес ма- 
териала, кг/м; И’ - работа сил трения, Дж; Их - скорость нарастания нагрузки, Н/с; 
К„з — Коэффициент взаимного перекрытия колеса и рельса; [, — интенсивность охлаждения 
колеса, Вт/м*; о — коэффициент теплоотдачи, Вт/(м".К); ^, — коэффициент теплопроводно- 
сти, Вт/(м.К); у — тепловой поток, Дж/с; И - износ весовой, кг/м?; А — амплитуда деформа- 
ции связей, м; й — шероховатость поверхности 2, м. 

Выполненная проверка по стандартным критериям гомохронности, Ньютона, Фруда, Био 
получила удовлетворительную сходимость полученных критериев подобия. Для подтверждения 
соотношения (4) составлен критерий экспериментальной триботермодинамики, позволяющий экс- 
периментально установить значение температуры вспышки Оьст активных микрообъемов: 

4 „О рой 
отрВАГу 
где а — показатель асимптоты градиента объемной температуры; т — масса контактная т,„, по- 

верхностная т„, либо объемная ти, кг (см. рис.1). 


Ко = = 14ет, 
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В целях прогнозирования образования ползунов на колесных парах, обусловленных тер- 
моповреждением активных микрообъемов колес при длительном вписывании в криволинейные 
участки пути или на сортировочных горках, были выполнены экспериментальные исследования по 
триботермодинамике. Эксплуатационные испытания (ВНИКТИ, 2007 г.) прохождения тягового ло- 
комотива в кривой малого радиуса по рельсам с установленной разомкнутой термопарой позво- 
лили установить поверхностную температуру контакта 580 и 637 °С [8]. Для подтверждения полу- 
ченных значений температуры и выдвинутой ранее гипотезы триботермодинамической модели 
фрикционного взаимодействия предлагается способ идентификации по трибоспектральным оцен- 
кам [9]. Исходными данными для идентификации выступали поверхностная или объемная темпе- 
ратура ®(#+х), регистрируемая термопарой, момент трения М(® и нормальная нагрузка МФ, сни- 
маемые с установки трения СМТ-1 при проведении эксперимента по схеме «тормозная колодка -— 
ролик». Внутри тормозной колодки была установлена хромель-алюминиевая термопара на задан- 
ном расстоянии от фактической площади касания. Начальный этап идентификации максимальных 
контактных температур и возможного термоповреждения заключается в использовании взаимной 
корреляционной функции (ВКФ) [10] как функции временного сдвига т 


В, (т)= | М(1)0(1+т)4. 
Анализ временных коррелограмм (рис.2) для всей реализации эксперимента позволил ус- 
тановить момент достижения корреляционной функции 0,8 И более, что может свидетельствовать 
о возможных термоповреждениях поверхностных слоев контакта. 
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Рис.2. Коррелограмма изменения взаимной корреляционной функции температуры и момента трения 


Определив временной диапазон анализа, далее выполняли анализ интегральных оценок 
степени диссипации по диапазонам частот, вычисленных по амплитудо-фазочастотным характе- 
ристикам [9], который показал, что потеря стабильности трибосистемы происходит с задержкой 
возрастания сил трения 5, причем градиенты изменения диссипации 3 и потери стабильности 4 
идентичны (рис.3). В частотном диапазоне 20-30 Гц фиксируются моменты формирования темпе- 
ратуры вспышки 11, а максимальное значение корреляционной функции отстает на величину пе- 
риода 12. 
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Рис.3. Соответствие интегральных оценок величины диссипации амплитудо-фазочастотных характеристик по диапазонам 
частот триботермодинамики фрикционного контакта термическому схватыванию [11]: 1 — максимум ВКФ температуры 
и момента трения; 2, 6, 9 — задержка распространения теплового потока; 3, 7, 10 — момент времени зарождения термопо- 
вреждения; 4, 8 — максимум потери стабильности движения; 5 — задержка распространения ВКФ; 11 — моменты времени 
теплоотдачи 


Анализ всех реализаций интегральных показателей диссипации трибосистемы показал, 
что наблюдаемое повышение величины диссипации трибосистемы с последующим развитием по- 
тери стабильности и повышения значения корреляционной функции характерно не всегда: после 
прекращения периода динамической нестабильности, возрастания диссипации трибосистемы и 
увеличения количества моментов потери стабильности, наблюдается возрастание корреляцион- 
ной функции, а задержка распространения теплового потока достигает одной секунды (рис.4). 

Оценкой формирования равновесной шероховатости и развития термоповреждений может 
выступать значение мощности трения: 

АР = тобьдот-т } 
где Р - сила трения, Н; Г - скорость относительного скольжения, м/с; /о(®, #) — интегральная 
оценка диссипации, вычисленная по амплитудо-фазочастотным характеристикам на задан- 
ных частотных интервалах со средней частотой ‹.,; и — частота вращения ролика; К - ра- 
диус ролика. 
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Рис.4. Соответствие интегральных оценок степени диссипации по диапазонам частот, показателя относительной неста- 
бильности трибосистемы и взаимной корреляционной функции температуры и момента трения атермическому схватыва- 
нию [11]: 1 — временной интервал разрушения равновесной шероховатости, увеличения инерционных и диссипативных 
сил и потери стабильности, образования локальных мостиков схватывания; 2 — формирование новой равновесной шерохо- 
ватости, повышение ВКФ температуры и момента трения; 3 — момент времени «зарождения» локальных мостиков схваты- 
вания; 4 — момент времени, при котором наблюдается «температура вспышки»; 5 — потеря стабильности по амплитуде 
или фазе; 6 — достижение ВКФ температуры и момента трения максимального значения, которое запаздывает на времен- 
ной интервал до 1 с 


На диаграммах зависимости мощности трения от объемной температуры фрикционного 
контакта на разных частотных диапазонах (рис.5) видно, что нормальное распределение мощно- 
сти трения характерно установлению равновесной шероховатости, отсутствию термоповрежде- 
ний, а экспоненциальное распределение — происходящим на фрикционном контакте термоповре- 
ждениям с выделением тепла. 
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Рис. 5. Оценка мощности сил трения по диапазонам частот как функция 
объемной температуры фрикционного контакта со средней частотой у: 
а-—5 Гц; б- 15 Гц; в- 25 Гц; г- 40 Гц 
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В целях определения максимальных значений «температуры вспышки» на фактической 
площади касания в соответствии с выдвинутой гипотезой экспериментальной триботермодинами- 
ки, были установлены в тормозной колодке уже две термопары на расстоянии 2 и 5 мм от факти- 
ческой площади касания с моделью колеса и зафиксированы две объемные температуры ©) и ©-. 
Результаты экспериментального определения поверхностной температуры на установке трения 
СМТ-1 представлены на рис.6. 
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Рис.6. Экспериментальная температура активных микрообъемов как функция сближения при скорости 
относительного скольжения: опыт 1 - 2,2%; опыт 2 - 1,5%; опыт 3 -— 0,5% скорости движения 


Для расчета температуры вспышки в соответствии с основными законами колебаний ме- 
ханических систем используем логарифмический декремент затухания [12]. Расчет выполняем 
следующим образом: 

1) темп затухания а, =©, /©, ‚ где ©, и ©) - значения температуры, измеренные термо- 


парами на разном удалении от фактической площади касания; 
2) зная логарифмический декремент ша; = .2, , ГДе 23, мкм — расстояние между установ- 
ленными термопарами, находим показатель асимптоты затухания а, = та. /2, ; 


3) согласно выдвинутой гипотезе экспериментальной триботермодинамики находим пока- 
затель асимптоты на номинальной площади касания <, =. -т,/т, , где ту — объемная масса; 
т», — поверхностная масса (см. рис.1); 

4) вычисляем логарифмический декремент затухания на номинальной площади касания 
Ша, = 9,2, , ГДЕ 2› - сближение поверхностей трения на номинальной площади касания; 


п! 


0222 __ 


5) исходя из выражения темпа затухания 4, =е”? =®, /©, находим поверхностную тем- 


пературу ®, = ©, /е“?. 

Зная величину сближения на фактической площади касания 21, аналогичным образом ус- 
танавливаем значение температуры вспышки ®©„„„, результаты расчетов которых представлены на 
рис.6б. Последующие вычисления температуры вспышки ®„„ по кривой опорной поверхности и 
предложенному алгоритму равенства произведений масс активных микрообъемов на показатель 
асимптоты подтверждают результаты эксплуатационных исследований поверхностной температу- 
ры, выполненных учеными ВНИКТИ [8]. С учетом предложенной методики триботермодинамики 
разработанная противоползунная конструкция ППС-12Д получила удовлетворительные эксплуа- 
тационные характеристики и рекомендуется для внедрения на всей сети ОАО РЖД. 

Заключение. В результате серии выполненных экспериментальных исследований триботермо- 
динамики по анализу амплитудо-фазочастотных характеристик, их интегральных показателей 
степени диссипации по диапазонам частот и степени нестабильности было выявлено следующее: 

— возрастание интегральных оценок степени диссипации более чем на величину полу- 
тора — трехкратного среднеквадратического отклонения относительно математического ожида- 
ния, вычисленного за продолжительный период времени, последующая регистрация потери ста- 
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бильности по амплитуде или фазе 4 (см. рис.3) и возрастание взаимной корреляционной функции 
6 при временном сдвиге, равном нулю, являются факторами, идентифицирующими возрастание 
поверхностной температуры активных микрообъемов соприкасающихся тел, образование местных 
мостиков сварки, флуктуации интегральных оценок диссипации 11, что характерно для термиче- 
ского схватывания поверхностей трения; 

— теплоотдача поверхностей трения характеризуется падением значений интегральных 
оценок степени диссипации и формированием новой равновесной шероховатости; 

— при атермическом схватывании (см. рис.4) было установлено, что локальные металличе- 
ские связи поверхностей трения обусловливают флуктуации интегральных оценок степени дисси- 
пации относительно их математического ожидания и сопровождаются потерей устойчивости по 
амплитуде или фазе амплитудо-фазочастотных характеристик 1, но не вызывают возрастание 
взаимной корреляционной функции температуры и момента трения; возрастание же последней 
характеристики 2 фиксируется после установления равновесной шероховатости и прекращении 
периода динамической нестабильности, т.е. задержка распространения теплового потока дости- 
гает одной и более секунды; 

— использование гипотезы экспериментальной триботермодинамики позволяет фиксиро- 
вать максимальные значения температур на фактической площади касания в текущем масштабе 
времени; 

— использование методик трибоспектральной идентификации и физико-математического 
моделирования позволяет установить идентификационные характеристики термоповреждений 
поверхностей трения подсистем «колесо — рельс», «колесо — тормозные колодки», что с исполь- 
зованием вибродиагностического комплекса по мониторингу термотрибодинамики позволит пре- 
дотвращать такие негативные явления, как образование ползунов на колесных парах, разрушение 
тормозных колодок и т.д. 
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ОУМАМТС МОМТТОВТМС ОР ЕКТСТТОМАЕ МОВТЕЕ $УЗТЕМ$ 
ТКЕТВОТНЕВМОБВУМАМТС$ 


А... ОГУАВКТМ 
(Во$оу Чае Тгапзрок Утмег$Ку) 


А ПуроПе${5 о! ехрептета! тБопегтодупат/с$ Ваа Бееп идде®еа апа ргоуеа. ТВе Вуроте55 регтй$ ю дЕ- 
пе Пар (етрегаёиге оп {пе гие сощасЕ агеа а$ а ргасйиса! ехатре. 

Кеуигога$: рпу/со-таетайса/ тоде/тд, ЕпБо5ресйа! а55еззтепЕ, Пазй етрегаиге, дпрртд, Гасе апа Бик 
[етрегаёиге$. 
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